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In trunking theory, peakedness is defined conventionally as the variance- 
to-mean ratio of a traffic load when carried on an infinite trunk group. 
For analysis of switching machine delays, it has proven useful to define 
a peakedness measure associated with the Carried Arrival Process (CAP), 
the stream of call arrivals carried on an incoming trunk group. The peaked- 
ness of the CAP is defined to be the conventional peakedness of a fictitious 
traffic-load process generated by associating with each carried arrival an 
independent exponentially distributed holding time with mean equal to 
the mean of calls actually carried on the trunk group. 

The problem considered is the effect of trunk group congestion on the 
peakedness of the CAP for traffic consisting of renewal inputs offered on 
a blocked-calls-cleared basis to a finite trunk group with exponential hold- 
ing times. The CAP is characterized as a semi-Markov process. This 
model leads to the determination of the peakedness of the CAP. Numerical 
results illustrate the reduction of peakedness, or smoothing, introduced 
by the congestion. 

I. INTRODUCTION 

This paper is concerned with characterizing the traffic offered to a 
switching machine, taking into account both the alternate routing 
that the traffic may have undergone and the smoothing of the traffic 
resulting from congestion on the trunk group incoming to the machine. 
In trunking theory, peakedness is defined conventionally as the 
variance-to-mean ratio of a traffic load carried on an infinite trunk 
group. It is well known that trunk group blocking of peaked traffic, 
such as overflow traffic, can be substantially larger than the blocking 
seen by Poisson traffic with the same intensity. Similarly, switching 
machine* delay and capacity can be quite sensitive to the peakedness 

* Throughout this paper, when we refer to a switching machine we mean the com- 
mon control devices in a switching machine. 
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of the incoming traffic. 1 To determine the peakedness of the traffic 
offered to a switching machine, we must take into account the smooth- 
ing effect of the incoming trunk groups. To this end, we consider the 
process of arrivals offered to a trunk group which are carried by that 
trunk group. We call this process the Carried Arrival Process, or CAP. 

To illustrate the CAP, consider the alternate routing network 
shown in Fig. 1. Here traffic overflowing trunk group A B is then 
offered to trunk group AC [Fig. 1(c)]. Those calls finding free circuits 
on AC then appear at node C as requests for service. The CAP is 
illustrated in Fig. 1(d). 

The basic model used in the analysis is shown in Fig. 2 where a re- 
newal process is offered to a group of N trunks on a blocked-calls- 
cleared (BCC) basis. The renewal input allows us to consider overflow 
traffic offered to an incoming trunk group. The holding times on the 
trunks are assumed to be independent, identically distributed, ex- 
ponential random variables with service rate /x- 

For analysis of machine performance it has proven useful to define 
a measure of peakedness for the CAP, z c , equal to the variance-to- 
mean ratio of the traffic load carried (number of trunks occupied) on 
an infinite trunk group to which the Carried Arrival Process has been 
offered. By definition, the holding times on the infinite trunk group 
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Fig. 1 — Carried Arrival Process. 
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Fig. 2 — Carried arrival model. 

are independent, identically distributed, exponential random variables 
(with service rate n) which are independent of the holding times on the 
incoming trunks. * 

The peakedness of the CAP should be distinguished from the 
variance-to-mean ratio of busy trunks on the incoming group, a 
quantity which is discussed, for example, in Section 8.4 of Reference 2 
for the case of Poisson input. This distinction can be made clear by 
considering an example of Poisson traffic of intensity X (calls/second) 
offered to N trunks. As X approaches zero, the two measures approach 
one. Clearly, as X gets large, the variance of busy circuits in the trunk 
group goes to zero and the mean goes to N, giving a variance-to-mean 
ratio of zero. On the other hand, as X gets large, the time differences 
between successive carried calls approach independent, exponential 
random variables with rate Nn (i.e., a Poisson stream) and the peaked- 
ness 2 C approaches unity. This is illustrated graphically in Fig. 3 
which plots z c and (v/m)u.s. (variance-to-mean ratio of busy servers 
on the N trunks) as a function of offered load for A r = 10. The example 
is a special case of the general results derived in this paper for arbitrary 
renewal input to the trunk group. 

By modeling the CAP as a semi-Markov process (SMP), it becomes 
possible to calculate peakedness z c as a function of the peakedness of 
the traffic offered to the incoming trunk group and of the congestion 
encountered on the group. The resulting z c may then be used in the 
determination of machine performance. 1 Numerical results illustrate 
the reduction in peakedness, or smoothing, introduced by the trunk 
group congestion. In the course of determining the peakedness, the 
transform of the distribution of the time between carried calls is 
derived. 



* That is, although carried arrivals are accepted simultaneously on the finite and 
infinite trunk groups, the departure times are different. 
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Fig. 3 — Distinction between peakedness and varianee-to-mean ratio of busy 
servers. 



II. OUTLINE OF RESULTS 

In this section, we give an informal overview of the results. In 
particular, we give the key equations that a user can employ to deter- 
mine the peakedness, z c , of a Carried Arrival Process. The equation 
numbers are the same as will be used in the derivation. Stationarity 
is assumed throughout. In all that follows, we assume unity holding 
time (or the time unit is the mean holding time). 

First of all, 

z c = M+-±, (11) 

where M + is the mean number of calls up on the infinite trunk just 
after the time that a call is accepted onto the finite (and infinite) 
trunk group. (1/0) is the mean number of carried calls (on both the 
finite and infinite trunk group).* 

M + is determined from the following equation : 

*+-r^-^[i-7(* + i)]- (32) 

where 0(s) is the Laplace-Stieltjes transform of the interarrival time 



'Note that, since we are assuming unity holding time, M + and (1//3) are in 
erlangs. 
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distribution. Mf is the mean number of calls up on the infinite trunk 
group immediately after arrival of a carried call, given j calls in progress 
on the finite trunk group immediately following the arrival of the 
carried call. P s (j = 0, 1, • • -, N) is the probability that an arrival 
finds j trunks busy on the finite access group, and 

pt - r^ (24)- 

is the probability of j calls up on the finite trunk group immediately 
following a carried call arrival. 

The only quantity left to be determined in (32) is M% which is 
calculated by solving the linear equations 



[M+ - 1]P+ = £ CiJItf", m = 1, 2, • • ., N, (33) 

Z=max(l,m— 1) 



where 



l^m^N, (34)+ 

p+ ( N \»-™+WN -m+l\ , 

( "- i - *(# + i) (35) 

A. simple method of solving (33) is discussed at the end of Section VI. 
In the course of deriving the expressions which ultimately determine 
the peakedness of the CAP, the Laplace-Stieltjes transform of the 
distribution of time between carried call arrivals is obtained. This is 
?iven by 

^-♦W-BT^Jfil*. (23) 

Note that the CAP is not completely characterized by (23), since it is 
not generally a renewal process. 
Examples are given in Section VII. 



* One method of computing P,- is via the equations given on p. 179 of Reference 3. 
Alternate methods which may avoid some of the numerical difficulties inherent in 
;his approach will be discussed in Appendix C. 

Alternate expressions for special cases, more suitable for computation, are given 
n Appendix C. 
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Fig. 4 — Peakedness of the Carried Arrival Process. 

III. CHARACTERIZATION OF THE CARRIED ARRIVAL PROCESS AS A SEMI- 
MARKOV PROCESS 

Consider a renewal process, with a nonlattice interarrival time dis- 
tribution F(t), as an input to N trunks with mutually independent 
exponential holding times each with mean unity (or the time unit is 
the mean holding time). ^(0 + ) is assumed to be zero. Blocked calls are 
cleared. Such a system is analyzed in Chapter 4 of Reference 3. 

Let {n, i = 1, 2, • • • } denote the sequence of times at which calls 
are accepted by the N servers (this is, of course, a subset of the times 
at which calls are offered). Let j(t) be the number of servers busy at 
time t. Then j(r+) is the number of servers busy just after the nth carry 
(nth carried call). Note that P{j(r+) = 0} = 0. It is clear that the 
j(t) process held fixed at j(r+) for r n < t < r n+ i is an SMR* The 
transition probabilities for the embedded Markov chain are derived 
as follows. 

Since there is a death process on the finite trunk group between 
carried calls, we have for m = 1, 2, ■ • ■, N and for I •» m — 1, • • •, 
N — 1 (in which case, the next arrival is the next carry) 

P{j(r+) =m|i(r+_ 1 ) = 1} 

= r ( l ^ e-^'Tl - e- l J- m+1 dF(t) 
J \m - 1/ 

=U- 1 )Tr: +1 ) ( - i) ^ +m - i) ' a) 

where <f>(s) is the Laplace-Stieltjes transform of Fit), 

4(8) = I" e-"dF(t)- (2) 

Jo 

* For an introduction to semi-Markov processes, see Reference 4, Chapter 5. 
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When i(T+_ x ) = JV, note that the next arrival need not be the next 
carry. With s.c. denoting service completion and F t (t) the e'th-fold 
convolution of F(t), we have 

P{j(r+) = m|j(r+_ 1 ) = N} 

fno s.c. before the first (i — 1) 
= .£ P j arrivals after r,*-!, (N - m + 1) s.c.'s j(r+_i) = N\ 
[before ith arrival after r n _i 



»=i Jo io \m — 1/ 

(."-.)TC"r +I )<-'>-<-+—> 



,-(m-l) (*-«)[" J _ e -((-a)-|iV-m+l 

XdF i ^ 1 (s)dF(t - s) 

N +1 _ 



0(iV) 



(3) 



Letting Fj M (0 be the conditional probability that a transition will 
take place within a time f, given that the process has just entered I 
and will next enter m, we have 



<tfW0 = 



(m - 1 ) ^"""'H 1 ~ ff-']*^ ,H "W(0 , 

P[j(r+) =m|j(r+_ 1 ) = Z} 

Z = m — 1, 

X[l-e-<'-«>]*-"+W,-- 1 (s)dF(<-a) 

^{i(r+) = m|i(r+_ 1 )=iV} 



<>-(m-l )(<-«) 



JV-1, 



J = iV. 



(4) 



Although the SMP characterization of the CAP is of general in- 
terest, it is particularly useful in determining the peakedness of the 
CAP as we shall see in Section IV. 

IV. PEAKEDNESS OF THE CARRIED ARRIVAL PROCESS 

Recall that the peakedness of a process is the variance-to-mean 
ratio of the number of calls up on an infinite trunk group when that 
process is offered to the infinite trunk group. To determine the peaked- 
ness of the CAP, consider the situation shown in Fig. 4. In this figure, 
each time a call is carried on the finite group, a call is put up on the 
infinite group with an exponentially distributed holding time with the 
same mean as on the iV-trunk group but independent of the AT-trunk 
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group holding time.* As before, j(t) is the number of busy servers on 
the N-trunk group at time t. Let k(t) be the number of busy servers 
on the infinite trunk group at time t. We use the following familiar 
result (given for renewal input in Reference 3 and for semi-Markov 
input in Reference 5, which is the form applicable to our problem). 
The following limits 

P» = lim P{k( T ~) =i\ (5) 



and 

exist and satisfy 



P; = \im P{k(t) = i\ (6) 



t-»00 



p* = ?Ltl (7) 

where /3 is the mean time between transitions of the SMP (i.e., mean 
time between carried calls). From (7) we obtain 

lim E{k*(t)} = I [ lim M(t-) + 1] - i lim M(r+) , (8) 

M(t) = E{k(t)}. (9) 

M+ = lim M(r+), (10) 



where we have defined 
Defining 



the peakedness of the CAP (denoted z c ) is given by 

2c = M+-|- (ll)t 

Note that 1/^ corresponds to the mean of the carried load (recall that 
we are assuming unity mean holding time). 

We are thus left with the problem of determining M + . This deter- 
mination will be in terms of the distribution of time between carried 
calls, to which the next section is devoted. 

V. DISTRIBUTION OF TIME BETWEEN CARRIED CALLS 

Consider an arrival at r n which finds a free circuit [i.e., j(r~) < iV]. 
Let F e (t) be the distribution of time until the next carry (carried call) ; 



* It is this independence which distinguishes the peakedness from the variance-to- 
mean ratio of busy trunks on the AT-trunk group (as discussed in Section I). 

f This is given in Reference 6 for the case of renewal input and weaker assumptions 
on service times. 
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i.e., 

F c {t) = P{ict^ Scarry at t„), 
where ict denotes inter-carry time. Denoting 

FcD\J(t+) <N1 = P{ict =g t\j(r+) < N] 
and 

F c {t) = P{ict^ t\j(r+) = N] 

and recognizing that 

^eD|i(r+) <N1 = F(t) 



yields 



where 



^.(»-to[i- t ^]+Aw[ i ^] i 



(12) 



(13) 

(14) 
(15) 



Pi = P\J trunks busy on the finite group just before 

a call arrival} (16) 

is the stationary call congestion probability given on p. 179 in Chapter 
4 of Reference 3. In particular, P N is the blocking probability. 
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Fig. 5— Peakedness of CAP (5 trunks). 
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Let <f> c {s) and c (s) be the Laplace-Stieltjes transforms of F c {t) and 
F c (t), respectively. Transforming (15) gives 

*.W^w[i-T^]+^)[l^]' W 

The function e («) can be obtained from the solution to the Type I 
counter problem given on p. 207 of Reference 3 with renewal input 
transform 0(s) : 

&(«) = [1 - *(«)] I™ e->»H(y)dm(y) , (18) 

Jo 

where, for our problem, 

H(t) = 1 - e~ Nl (19) 

and 

m(t) - E\ number of arrivals in (0, t)\. (20) 
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Fig. 7— Peakedness of CAP (20 trunks). 
From (18) we have 

where we have used 



4.0 



/; 



e-"dm(t) = r ^ n 
1 - 0(a) 



(21) 



(22) 



Combining (17) and (21) gives the transform of the intercarry time 
distribution : 

We are now in position to determine M + , defined by (10), and sub- 
sequently the peakedness of the CAP. 
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Fig. 8— Peakedness of CAP (50 trunks). 
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VI. DETERMINATION OF M + 

In order to determine the peakedness of the CAP [eq. (11)], we 
need to evaluate M+ defined in (10). This will be done by characterizing 
the conditional mean of trunks up on the infinite trunk group, given j 
calls up on the finite trunk group. 

Recall that we are considering an arrival at r„ that finds a free 
circuit on the finite trunk group (i.e., j(t~) < N). The state dis- 
tribution on the finite trunk group at t+ is thus given by 

Pf- = Pr{ j{it) - i) - f^ . !****> (24) 

Pf - , 

where the P/s are the call congestion probabilities defined in (16). 
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Fig. 9— Peakedness of CAP (100 trunks). 
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We define 



and 



P$ N = Pr{j(T+) <N\, 

M~ = E\k(r-)} = E{k(r+)\ - 1 = M+ - 1, 
M+ = E{k(r+)\j(r+) =1], 

Mt N = E{k(r+)\j(r+) <N\, 



(25) 
(26) 
(27) 

(28) 

where k corresponds to the infinite trunk group and j corresponds to 
the finite trunk group. In terms of these quantities, we have 

M+ = M+ N P+ N + M%Pt. (29) 

Recall that, if j(r+) < N, the distribution of time until the next 
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0.1 0.2 

Fig. 10 — Peakedness of CAP versus blocking (5 trunks). 

carry is F(t) and if j(r+) = N the distribution of time until the 
next carry is F c (t). Using this together with some conditioning argu- 
ments, the following relationship is obtained (see Appendix A) : 

M~ = M+-l= P% N Mt N <t>(l) + PtMUM). (30) 

From (29) and (30) we obtain 

+ *&?[_ 1-0(1) J 



M+ = 



1 - 0(D 



(31) 



Use of (21) simplifies (31) to 



M+ =T 



1 



M+PU(N + 1) 
1 - 0(1) 1 - <f>(N + 1) 



(32) 



It should be noted that the first term in (32) corresponds to the 
value M+ would assume if the renewal input process was offered 
directly to the infinite trunk group. The second term corresponds to 
the reduction in M+ as a result of blocking on the finite trunk group. 
We are now left with the problem of determining M %. 

It is shown in Appendix B that M+ for in = 1, 2, ■■■, N satisfies 
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Fig. 11— Peakedness of CAP versus blocking (10 trunks). 

the following set of equations : 

[M+ - 1]P+ = £ C lm Mf, m « 1, 2, ■ • -, N, (33) 

l =max(l,jn — 1) 

where for m — 1 fS I < N — 1, 1 ^ »i g N, 

*■"•+» // - m + 1 

Further, for I = N we have 

JV \^-^i+i /AT - m+ 1 



c - -*■(»! OX (' ; ,+1 )(-D'^+»)- (34) 

P Hm-l) ,5 ( , )(-DV(* + i) 



VVm — 



(35) 



1 - *(tf + 1) 

It should be noted that the above set of equations is in a form which 
is amenable to solution for the desired quantity M%. Written in 
matrix-vector form, the matrix in question is triangular with additional 
entries below the diagonal. Transforming the matrix associated with 
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Fig. 12— Peakedness of CAP versus blocking (20 trunks). 

(33) into triangular form* leads directly to the quantity of interest 
M% for use in (32), which is subsequently used to determine the 
peakedness of the CAP [eq. (11)]. 

VII. EXAMPLES 

We ran some examples with a 2-moment match 1 interrupted Poisson 
process (Reference 7) as the renewal input to the finite trunk group 
(the computational aspects are discussed in Appendix C). Figures 5 
to 9 show z c , the peakedness of the CAP, as a function of p, the offered 
load per trunk for N = 5, 10, 20, 50, and 100 trunks, respectively. 

* Details are in Reference 9. 

+ The blocking experienced by an overflow process is less than the blocking seen 
by a 2-moment match interrupted Poisson process and more than that seen by the 
3-moment match process (all with the same mean and peakedness). 
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Fig. 13 — Peakedness of CAP versus blocking (50 trunks). 

On each figure are plots for offered z's of 1, 2, 3, and 4. It is seen that, 
for fixed N, as z increases, the smoothing effect (reduction of peaked- 
ness) becomes appreciable at lower p's. This is because of blocking 
remaining negligible for larger p's as z is decreased. If we fix z, we see 
that the smoothing effect becomes appreciable at lower p's as N is 
decreased. This is again explainable from the point of view of blocking, 
i.e., blocking is larger on the less efficient small trunk groups. 

Since blocking is an important parameter, Figs. 10 to 14 show the 
peakedness of the CAP versus the blocking for the same cases as 
shown in Figs. 5 to 9. Note that, for final trunk groups which are 
normally operated with blockings of 0.01, the smoothing effect is very 
small, while for high-usage trunk groups which may reach blockings 
of a few tenths, the smoothing is substantial. Also, note that z c , in all 
the cases, approaches unity as the load (and blocking) increases which 
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Fig. 14 — Peakedness of CAP versus blocking (100 trunks). 

is consistent with the explanation in Section I as to how peakedness 
differs from the variance-to-mean ratio of busy servers (which ap- 
proaches zero). 

It was observed in Reference 10 that, when z > 1, the blocking 
probability is bounded away from zero no matter how small the input 
mean is. This is evident in Fig. 10. 

VIII. CONCLUSION 

We have shown how to determine the peakedness of a CAP. The 
use of mean and peakedness to characterize a CAP is attractive from 
the point of view of simplicity and is consistent with the use of the 
equivalent random method (Reference 2) in trunking analyses. To 
approximately calculate the delays at a switching machine, we could 
replace the CAP (or, more usually, a superposition of CAP's) with an 
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interrupted Poisson process with the same mean and variance and 
proceed as in Reference 1 . This is being investigated. 

We mention, in passing, that we also tried out a renewal approxi- 
mation for z c . That is, although the CAP is a semi-Markov process 
and not generally a renewal process, we tried to approximate z c with 

1 1 



l-0 e (l) £' 

which is the peakedness that the CAP would have if it were a renewal 
process. This approximation did not compare well enough with the 
true z c to recommend its use. That is, although one might use a re- 
newal approximation to a superposition of CAP's after the peakedness 
is determined (sec the last paragraph), we do not recommend using a 
renewal assumption to determine the peakedness. 

In the course of determining the CAP peakedness, we have more 
fully characterized the CAP as a semi-Markov process. Any queuing 
results available for semi-Markov inputs could be used with the CAP 
semi-Markov characterization given in Section III. 
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APPENDIX A 

Derivation of Equation (30) 

Consider the system in equilibrium with times of carried calls {r n j. 
From (26) we have 

M~ = E{k(r-)\, (36) 

which can be expanded as 

M- - Eik(r-)\j(r+. 1 ) = N}P{j(r+. 1 ) = N] 

+E{k(r-)\j(r+_ 1 ) < N\P{j(t+-. 1 ) < N\. (37) 

This can be written as 

M- = P$±± iP{k(r-) = t|t(i* = n, HriLt) = N\ 
XP{fc(r+_ 1 ) -nljfoki) =N] 
+P$n t t iP{k(r-) = i\k(r+^) - n, j{rt-x) < N\ 

n=0 t'«4 

XPik(T+- 1 )=n\j(T+- 1 )<N} . (38) 
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Writing 

P{k(r~) = z|fc(r+_ 1 ) = n, j(r+.,) = N} 

- f*P{*(K) = i\k(rt-J - ~, i(rti) = iV, 

r„ - r n _ x = *}dP e (0 

= r (?) <r"(l - e-'^-HFcit) (39) 

and 

P{/b(r-) = i|fc(r+_,) = W, j(lt-l) < N] 

= f° (?) e _it (l - (T^-'dFit) (40) 
and observing that 

X) i ( ?) «-*'(! ~ e-') w_ * = we- 4 (41) 

leads to the desired result 

M- = PfcrJIti^a) + PtMUM. (42) 

APPENDIX B 

Derivation of Equations (33) through (35) for M+ 

Consider the events 

Ai m = {j(-r^-i) = I, (I — m + 1) s.c.'s before next 

arrival after t„_i} (43) 

B km = { j(r£.i) = iV, no s.c. before (k — l)st arrival after r B -i, 

(JV — m H- 1) s.c.'s before A;th arrival after t„_i}, (44) 

where s.c. denotes service completion on the finite trunk group. From 
these definitions we obtain 

P{j(r+) = m] = L 1 P{A lm ] + t P{B km \. (45) 

The conditional mean of interest is thus given by 
M~ = Eik(r-)\j(r+)=m} 

N £ E{k(r-)\A lm }P{A lm \ + t E{k(r-)\B km }P{B km \ 

_ l=rn-l k=l (46) • 

Defining the events 

Ai={j{rt- l ) = l) (47) 



* Whenm = 1, P{A 0m ) = 0, because P{ j'(t+_j) =0} = 0. In that case, sum from 
I = ltol = N -1. 
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and 

Au m = [ (I — m + 1) s.c./s before next arrival after r„_i} (48) 
gives 

E\k(r-)\A lm \P{A lm } 

= t iPik(r-) = i, A 2 i m \A l }P+ (49) 

= P,+ t P{*(rt,) = r|A,j £ tP{fc(r a -) = t, 

r-=0 ,=1 

^2 im |A ; , fc(T+.,) = r\. (50) 
Letting Z < N and using 

**{*(*£") = •', A Um \Ai, fc(r+_i) = r) = /"" ( r ) e -«'(l - e -') r -'' 

X (m - 1 ) e_(m_1) '(l - e->y->»+HF{t) (51) 
and 

.£ i ( J j e~»(l - e-<)'-« = re-' (52) 

yields 

E{k{r-)\A lm }P\A lm \ = P+M+ J* ( m l _ \ 

Xe-»'[1 - e-tJ-n+WFit). (53) 
From (1) we have 

E{k(r-)\A lm )P\A lm \ =P+M?( l _\ 

X £ (' * +I )(-D^ + ») (54) 

= Ci m ilf f +, (55) 

where we are using (54) to define C [m in (55). This yields (34). 

It now remains to show the desired relationship for C Nm . We con- 
sider the second summation of (46) 

E{k(r-)\B km \P\B km ) = tiPik(r-) = i,B l ,B 2hm } ) (56) 
where 

Bi= (i(r+_ 1 ) =N] 
B 2km = {no s.c. before (k - l)st arrival after r n _i, (N - m + 1) 

s.c.'s before ifcth arrival after t„_i}. (57) 
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Now 
E{k{r-)\B km \P{B km \ 

= Z t *P{*(t« ) = h B ikm \Bi, fc(T+.j) = n\ 

XP{fc(r+_ 1 ) = n\B 1 }Pt (58) 
n i" 00 f* n / n — 1\ 

X(l - rO'-V" ( ^ \«-<»-i><«— ) 

x (1 _ e -(^.))Jv-» +HF k -x(s)dF(t - s). (59) 
Upon performing the summation over i, (59) simplifies to 
-E{fc(r„-)|.B* m }P{£* ra } 

X N 'f: +1 ( N ~ m + 1 ) (-iye-^-'HF k - l (s)dF(t - s) 
v to \ D / 

/ jy \ i\r-m+i /JV-m+l\ 

= W(/l) & ( - X) '( t ) 

X**" 1 ^ + l)*(u + m). (60) 
Summing over fc gives 

t E{k{r-)\B km \P{B km ) 



k=\ 



n( m !,) t'H)f-; +1 )^^» 



M+ "\w-l/ ,=o \ 5 / (61) 

- M * 1 - *(JV + 1) 

which is the desired result for (35). 

appendix c 

Computational Considerations 

In this appendix we briefly* discuss some computational problems 
experienced in the numerical solution of the carried process problem 
and point out some possible approaches to circumvent them. f In 



* A more detailed description of our computational experience is in Reference 9. 
+ We first discuss the approaches and then the numerical experience we have nad. 
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particular, we are concerned with the computation of the state prob- 
abilities, P h j = 0, ■■■, N, defined in (16) and the computation of 
the quantities Ci m defined in (34) and (35). 

We also specialize and subsequently simplify the above results for 
the case where the interarrival time distribution is a sum of ex- 
ponentials. In particular, we consider the interrupted Poisson process 7 
as the renewal input. This specialization appears to have eliminated the 
numerical problems associated with the required calculations for 
large trunk groups. 

The problem of determining the call congestion state probabilities 
for a renewal input to a BCC system with N independent exponential 
servers is considered on p. 179 of Reference 3. The results are as 
follows : Let 

^-( llrtJ + l) ) " i-%-.X-l. (02) 
with Co = 1. Then B r) the rth binomial moment of the P/s, is given by 

The B r satisfies the backward recursion 

*-(4#S fl )*~+(?)*'' w 

1 1 " /N\ 1 

RZ = p~ = S U- ) r~ 65 

(Reference 11, p. 93). The P/s are given by 

p i= r t i-iy-t^Br- (66) 

The computation of the Cj coefficients and the binomial moments 
B r is fairly straightforward and does not pose much of a numerical 
problem. It is the computation of the P/s, using (66), that is sensitive 
to numerical errors. The alternating sign in (66), together with the facts 

that ( J B r can be quite large and the summation in (66) is between 

zero and unity, lead to a numerical problem.* 

* Actually, for an N = 10 case seven significant decimal digits were lost in one 
subtraction and the resultant probability was computed to be zero. This plays havoc 
with the solution to (33). Details are in Reference 9. 



with 
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An alternate approach* to the computation of the state prob- 
abilities is by way of the equations 

P k = £ PikP,; (67) 

where the transition probabilities p jk represents the probability of 
going from state j prior to one arrival to state k prior to the next 
arrival. The p/*'s are given by 

Vjk . { * + * W" e-*'(l - e-<y+i-WF{t) , j <N, (68) 

and 

Ptf.fc = PN-i,k- (69) 

From (67) and (68), we obtain the backward relation 

F "-m-m& Pi>p '- (70) 

Expanding (68) gives 

*»=( j tT£( j+ l- k )^ +k) - (71) 

The computational procedure is outlined as follows: Compute P N 
from (65) f and use (70) to compute P t for j < N. It should be noted 
that, although the terms in (71) alternate in sign, <f> never exceeds 
unity and is monotonically decreasing. Also note the relation between 
p jk given by (71) and C im given by (34) and (35). At this point, the 
accuracy in computing the P/s should be comparable to the accuracy 
in computing the Ci m 's. 

For the case where F(t) is the sum of exponentials (e.g., interrupted 
Poisson process) we can further simplify (and more accurately com- 
pute) the P/s and C tm . We go to the integrals from which the sums 
(with alternating signs) appearing in (34), (35), and (71) were derived. 
Note that we have 

( 1 1 )-T( , — +1 )(-i)^+«) 

= f ( m i i ) e ~ mt ^ - ^y-^HF{t) (72) 

* Motivation for investigating this approach stems from remarks made by P. J. 
Burke. (In recent unpublished work, Burke showed a more accurate approach for 
the case where the renewal interarrival time distribution is a sum of exponentials.) 

+ Note that each term in the sum of (65) is positive. 



PEAKEDNESS OF CARRIED ARRIVAL PROCESS 1641 

[see (53) and (54)]. Let 

F{t) = E fc.(l -«-*); (73) 



i=i 



then the integral in (72) can be identified as a beta function.* Re- 
peated integration by parts in (72) gives 



/o" ( m - 1 ) e_m ' (1 - r ') ! " n+W W " /£/■*<** r,), (74) 
where 

/ Z - (Z - m) \ , , 

"A i + r,-g-») J for *>»-!. (75) 

For I = m — 1, we obtain from (72) 

/—!..(*« r«) -~|~ (76) 

Note that / ifn can be computed recursively from 

/l+i.«(fc<, r.) = ( ; + r , + 2 )/*..» (k, r.) (77) 

with initialization from (76). 

The direct calculation of the integral has thus led to a computa- 
tionally tractable method of computing the Ci m 's and the P/s. C !m is 
calculated from 

C lm = P+ ^ ZfUK n)\ , m - 1 ^ I £ N - 1 , (78)t 

[see (53)] and 

CArm = 1 - </>(# + 1) ( £/*■-(** r <) J (79) f 

[see (35) and (72)]. 

The P/s are calculated from (70) where (68) is simplified as above 
(using integration by parts) and then used to compute the transition 
probabilities. Note that for an interrupted Poisson process, 7 S = 2, 
and hi, fc 2 , r h and r 2 are given in Reference 7 in terms of the switch 

* Identification made by D. L. Jagerman. 
Note that this procedure does not involve the calculation of either binomial 
moments or binomial coefficients. 
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parameters. An alternate method of computing the P/s for the 
interrupted Poisson process is via the use of birth and death equations 
and conditioning the results on the switch being closed. A computer 
program for doing this was available (Reference 8). 

At this point, it is of interest to discuss our computational experi- 
ence using some of the aforementioned procedures. The first method 
considered was to calculate the P/s by first obtaining the binomial 
moments [eqs. (62) to (66)] and to compute the Ci m 's from (34) and 
(35). Using single precision arithmetic the procedure worked for N = 2, 
but failed for N = 10. The problem was traced to inaccurate calcula- 
tion of the probabilities from binomial moments. Double precision 
arithmetic extended the range of N (worked for N = 10). The method 
failed at N = 20. The failure was traced to the same cause as above. 
At this point we used the birth and death equation approach to calcu- 
late the P/s, 8 which assumes an interrupted Poisson input. This 
extended the range to N = 20. For N = 30 we ran into problems 
computing the C im 's from (34) and (35).* Modification of the C tm 
computation using (78) and (79) significantly extended the useful 
range on N. The results presented in Section VII were computed using 
this method of calculation. 
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